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Abstract 



q 

O . This paper studies tliree types of functions arising separately in the anal- 

ysis of algorithms that we analyze exactly using similar Mellin transform 
techniques. 

> 

' The first is the solution to a Multidimensional Divide-and-Conquer (MDC) 

, recurrence that arises when solving problems on points in d-dimensional 

space. 



The second involves weighted digital sums. Write n in its binary represen- 
' tation n = {bibi^i •••6160)2 and set SMin) — J2t=o^^^^t'^* ■ analyze 

O . the average TSuin) = ^ I]j<n Smij)- 

The third is a different variant of weighted digital sums. Write n as n = 
2n 2""^ + ■ ■ ■ + 2*" with ii > Z2 > • • • > ifc > and set WM{n) = 
X ■ t^'TK We analyze the average TWM{n) = i E,<„ WmU)- 



We show that both the MDC functions and TSuin) (with d ^ M + I) 
have solutions of the form 

d-2 

XdTiXg'''^^ n + ^ (nlg™ri)yld,m(lgn) + Cd, 

where A^, Cd are constants and Ad^miuy^ are periodic functions with period 
one (given by absolutely convergent Fourier series). We also show that 
TWM{n) has a solution of the form 

Af-l 

nGM (Ig n) + dM Ig'^n+Y, (ig'' GM.d (Ig n), 

d=0 
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where cIm is a constant, Gm{u) and G'A/,d(w)'s are again periodic functions 
with period one (given by absohitely convergent Fourier series). 

1 Introduction 

In this paper we use Mellin Transform techniques to analyze three types of 
functions arising separately in the analysis of algorithms: Multidimensional 
Divide-and-Conquer and two different types of weighted digital sums. 

(A) Multidimensional Divide-and-Conquer: 

The Multidimensional Divide-and-Conquer (MDC) recurrence first appeared in 
the description of the running time of algorithms for finding maximal points in 
multidimensional space. Previous analyses by Monier [18] gave only first order 
asymptotic, showing that the running time for the d-dimensional version of the 
problem is (Ig n = log2 n) 

Td{n) = XdU Ig'^-^ n + o{n Ig-^'^ n) 

for some constant A^. We will extend the Mellin Transform techniques for 
solving divide-and-conquer problems originally developed in |10] (see [13] for a 
review of more recent innovations) to derive exact solutions, which will be in 
the form of 

Td{n) = XdU Ig^-i n + ^ (n Ig'" n) ^rf,^(lg n) + c^, (1) 

m=0 

where A^, q are constants and m(n)'s are periodic functions with period one 
given by absolutely convergent Fourier series. 

(B) Weighted Digital Sums of the First Type: 

The second type of function we study is a generalization of weighted digital 
sums (WDS). Start by representing integer n in binary as n = {bibi^i ■ ■ ■ 6160)2- 
Define 

i 

Si{n) ■.= Y,m\ 

t=0 

i.e., weight the t^^ digit by its location in the representation. One can also 
view this as 2 times X](=o '^^ich is analogous to the derivative of the 

binary representation of re. This sum arises naturally in the analysis of binomial 
queues where Brown [6] gave upper and lower bounds 

[re Ig n — 2re] < Si (n) < [re Ig rej . 

Generalizing 5*1 (n) allows the "weights" to be any polynomial of t. Set 
5o(re) := n and VM > 1, define 

i 

SM{n):=^t^bt2\ (2) 
t=o 
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Figure 1: The graphs of (S'A/(n) — nlg*^n)/(nlg^^~^ n) for M = 1 (left) and 
M = 2 (right) plotted against Ign. Although the functions appear periodic 
they possess "large" fluctuations. These make direct analysis of SM{n) difficult, 
suggesting the analysis of its smoothed average instead. 



where t^^ := t{t + l)(t + 2) ■ ■ ■ {t + M - 1) is the M*^ rising factorial of t. 
The function SM{n) is not smooth (see Figure [T]). We will instead analyze its 
average 



TSM{n) 



1 



n 



j<n 



SmU)- 



(3) 



We will show that, surprisingly, TS'jvf(?7-) has an exact formula, which is in 
exactly the same form ([1]) derived above for the MDC problem (with d = M + 1 
and different constants). 

(C) Weighted Digital Sums of the Second Type: 

The third type of function we study is another WDS variant. Its simplest form 
arises when analyzing the worst-case running time of bottom-up mergesort. 
Assum^ that the worst-case running time to merge two sorted lists of sizes ni 
and 77-2 into one sorted list is rii + n2- 

Define C^(n) to be the worst-case running time of bottom-up mergesort 
with n elements. Bottom-up mergesort essentially splits a list of n items into 
two sublists, sorts each recursively, and then merges them back together. If n 
is a power of 2, then it splits the list into two even parts. If n is not a power 
of 2 though, i.e., n = 2^ + j with 1 < i < 2^^ — 1, then the algorithm splits the 



items into one list of size 2^ , 
satisfies the recurrences: 



one list of size j. Thus it is known that Cw{n) 



CM'') 



A;2* 



C^(2'=) + CUj) + (2' + j), for 1 < j < 2^ - 1. 



^The actual worst-case time is ni +n2 — 1. But, any mergesort uses exactly n — 1 merges, 
so the running time derived with cost ni + n2 is exactly n — 1 more than the real worst-case 
running time. 
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Panny and Prodinger [19] derived an exact solution for C^(n) containing a 
term G(logn), where G{x), defined by a Fourier series, is periodic with period 
one. However, the Fourier series is only Cesaro summable. Furthermore, G{x) 
is discontinuous at all dyadic points (points of the form x = n/2™', where n is 
integer, m is non- negative integer), which are exactly the points of interest. 

In this paper we will decompose Cw{n) into two different types of WDS 
and analyze the (smoothed) average of each part. We will then generalize the 
functions found and analyze the generalizations. 

Starting with the binary representation of n, ignore the bits and write 
n as the sum of descending powers of 2, i.e. n = 2*^ + 2*^ + . . . + 2'* with 
zi > ^2 > • • • > ^fe > 0. Iterating the above recurrence for Cw{n) gives 

k k k 

Cw{n) = ^t2** + Yl - 2*^- = Si{n) + ^ t2** - 2'\ 
t=i t=i t=i 

where Si{n) is the WDS of the first type defined previously. This motivates 
the introduction of another variant of WDS: 

k 

lyi(n) := ^t2**. (4) 

t=i 

As with the WDS of the first type, Wi{n) is not smooth enough to be analyzed 
directly (see Figure [2]) , so we instead study its average 

TT^i(n) := i^H^i(j). (5) 

j<n 

Similar to the WDS of the first type, this problem may be generalized by 
weighting the powers of 2 with polynomial weig htH, i.e. by defining Wo{n) := n 
and, VM > 1, 

k 

WM{n):=Yt^'2'^ (6) 

t=i 

and then introducing the average functions 

TWM{n):=^Y^MU)- (7) 

j<n 

We will show that TWMin) has an exact closed- form formula, which is in 
the form of 

M-l 

TWMin) = nGMilgn) + du \g^' n+Y, ^) GA/,d(lgn), (8) 

^ For WDS of the first type we use weights of the form t''' ; for WDS of the second type 
the weights are of the form t^' . The difference is due to ease of analysis. Both types of span 
the space of polynomials and hence our study allows any polynomial weights. 
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Figure 2: WM(n)/n plotted against Ign for M = 1 (left) and M = 2 (right). 
Though these functions also appear periodic, they are not smooth and are with 
"large" fluctuations. These make direct analysis of WMin) hard. 

where cIm is a constant, Gm('u) and G'Af,d(tt)'s are periodic functions with period 
one given by absolutely convergent Fourier series. 

Our approach to solving all three problems will be similar. We first use the 
Mellin-Perron formula and problem-specific facts to reduce the analysis to the 
calculation of an integral of the form 



3+ioo 



K{s)ds 



3— ioo 



for some problem specific kernel K(s). We then identify the singularities and 
residues of K{s) and use the Cauchy residue theorem in the limit to evaluate 
the integral. 



2 Background 



2.1 The Mellin-Perron Promula 



The main tools used in this paper are Dirichlet generating functions and the 
Mellin-Perron formula. For more background see, [20^ pp. 13-23], [9j and [12^ 
pp.762- 767]. 

Theorem 1 (The Mellin-Perron formula). Let {Xj}, j = 1,2, . . . be a sequence 
and c> lie in the half-plane of absolute convergence of Yl'jLi ■ Then for 
any m > 1, 



1 y-A.ri-iV"=^r"ff ^ 



j<n 



n 



27ri 



n^ds 



^ j s{s + l){s + 2)---{s + m) 



(9) 
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6 7 8 I) 10 11 12 13 14 6 7 8 9 10 11 12 13 14 



(e) TWi{n)/nvs. Ign (f) TW2{n)/nvs. Ign 

Figure 3: These figures illustrate the periodic nature of the second order asymp- 
totics of (when k = 3,4), TSuin) (when M = 1,2) and the first order 
asymptotics of TWMin) (when M = 1,2). 
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In particular, when m = 1 and m = 2, 



1 , . . 1 r^'°° I ^ X-i \ n'ds 



n ^ — ' 





j-c+ioo 


2m . 


1 c—ioo 




/•c+ioo 


27d . 


J c—ioo 



We will analyze the MDC functions and the two types of WDS by rewriting 
them to summations in the form of the left hand side of (llOp . WDS of the second 
type will also need a summation as in the left hand side of (|lip . The Mellin- 
Perron formula will then enable us to evaluate the associated line integrals 
instead. The line integrals will be evaluated exactly via the Cauchy residue 
theorem by considering integrations over some special contours. 

In the right hand side of ([1]), X^j^i Xjj~^, the Dirichlet generating function 
(DGF) of {Aj}, is the only factor depending upon {Aj}. The Cauchy residue 
theorem relates the value of the line integral to the residues at the poles of 
the kernel in the line integral, thus understanding the locations and associ- 
ated residues of the DGF's singularities will be essential to evaluating the line 
integral. 

Define the backward difference function VA by V A[j) = A{j) — A{j — 1) 
for any function A. The following lemma will be needed later in the analysis of 
WDS. 

Lemma 1. Let A he a function with A{{)) = and 



n 



Then 



1 /•^+*°° l^VA{j)\ n'ds 



where c > lies in the half-plane of absolute convergence of Y2'^^'\/A{j)j~'^ . 
Proof Note that TA{n) = ^ J2j<ni^ - i)V^(j) and then apply 1^0^. □ 

A similar lemma, previously proven by Flajolet and Golin [IQJ, will be 
needed to analyze the MDC functions. For any sequence {s„}, define its double 
difference sequence {AVs„} by AVs^ := -Sn+i — 2s„ + s„_i for all n. 

Lemma 2. Consider the recurrence 

fn = /[n/2j + fln/2] + (12) 
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with boundary conditions eo = ei = and /i = 0. Then 



and 

where c lies in the half-plane of absolute convergence of Y^J^^ AS/ejj~^ . 

2.2 Useful Facts Involving the Riemann-Zeta Function 

The Riemann-Zeta function is defined by C{s) := X]i>o^ when Re{s) > 1. 
Since it will appear in the integral kernels in the analyses of the WDS, we list 
some basic facts concerning the Riemann-Zeta function \23\ [24] that we will 
need. 

First, C{s) can be analytically continued to be analytic in the whole complex 
plane with the exception of a simple pole at s = 1 with residue 1. 

Next, in [9], Flajolet et. al. proved the identity 

/ Cis)—, r = 0. (14) 

By mimicking their proof, we prove the similar formula (for completeness the 
proof is provided in Appendix [B]) : 

1 /-5/4+*~ n'ds 
^ C(s)^— T^— =0. (15) 
27rz y_5/4_joo s{s+l){s + 2) 

When integrating C(s) the following asymptotic bounds [23] will be useful: 

Lemma 3. If s = a + it, where cr, t € M, the Riemann-Zeta function satisfies 
the bound 

c(s) = o(|trMiog|t|) (16) 

where 

\ — CF for fj < 



\ forO<.<i 

l-afor\<a<l 
/or 1 < a. 
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2.3 Useful Formulae Involving Some DGFs 

To understand the locations and associated residues of the integral kernels, we 
will need closed-form formulae of their associated DGFs. We start with some 
basic definitions. 

Definition 1. Express n = (bibi^i ■ ■ ■ 6160)2 in its binary representation. 
Set v{n) := X]t=o ^* number of "l"s in the binary representation of n 

and V2{n) to be the number of trailing "0"s in the binary representation of n. 
For example, 

ifn = U = (101100)2 then v{n) = 3 and V2{n) = 2; 
i/n = 33 = (100001)2 then v{n) = 2 and t;2(n) = 0. 

We can now introduce two useful DGFs. 
Definition 2. VAf > 0, denote the DGFs of v{n)^'^ and {v{n) + V2{n))^ by 



The analysis of these DGFs will require the following facts. 
Lemma 4. Let n be a positive integer. Then 

1. v{2n) = v[n) and v{2n + 1) = v{n) + 1; 

2. V2i2n) = V2{n) + 1; 

3. if n is odd, V2{n) = 0; 

4. v{n) — v{n — 1) = 1 — V2{n). 

Proof. If n = (6i6i_i ••• 6160)2 then 2n = (6j6j_i ••• 6160, 0)2 and 2n + 1 = 
(6j6j_i • • • 6160, 1)2, so v{2n) = v{n) and v{2n + 1) = v{n) + 1. 

If the binary representation of n has t trailing "0"s, then the binary repre- 
sentation of 2n will have t -|- 1 trailing "0"s. This proves V2{2n) = V2{n) -|- 1. 

If n is odd, the rightmost digit of the binary representation of n must be 
1, i.e. there is no trailing "0" in the representation. Hence V2{n) = for odd 
integer n. 

Slightly rewriting n as n = (6j6j_i6t+i, 1, 0, • • • 0, 0)2 where t = V2{n) shows 
that n — 1 = (6j6i_i6t+i, 0, 1, 1 • • • 1, 1)2. Therefore Vt>(n) = v{n) — v{n — 1) = 
1 - V2{n). □ 

From Lemma m it is straightforward to prove the following lemma, which 
includes formulae expressing some special DGFs in terms of C,{s) and Vm{s). 
For completeness, we provide its proof in Appendix [Ul 
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Lemma 5. For M > 1, 



odd j 

The following DGFs have closed-form formulae in terms of C,{s): 

t'-f~,as). (19) 

.7 = 1 



2.4 Absolute Convergence of Fourier Series 

In all three problems, evaluating the line integrals of the kernels will reduce to 
summations of residues at poles regularly spaced along a vertical line. These 
summations will best be expressed as Fourier series. To be useful, we will need 
to show that these Fourier series converge absolutely. Our major tools will be 
the following two lemmas. 

Lemma 6. Let e > 0, o"o,fo G M, > 1 + e o,i^d f be a complex function. If 

1. f is analytic in X = {s = a + it : a > uq — e,\t\ > to ~ e} o-nd 

2. 3A,B>0 suchthatya + iteX, \f {a + it)\ = 0{\t\^ \og^ \t\), 

then, for every fixed integer q > 0, 

Va > (70, V|t| > to, + it)\ = 0{\t\^ log^ 



Proof. From the Cauchy integral formula, for all s = a + it with cr > ctq and 
\t\ > to, 



2m J iz- 



where C = {z :\z — s\ = e}. Hence 



{z - 



dz 



< ^x{2ne)xO{\t\^log^\t\)x-^ 

= 0(|t|-^log^|t|) 



for fixed q and e. 



□ 
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Before stating the next lemma, we clarify that the statement has a 

pole of order at most at s = sq" ) allows the possibility that h{s) is analytic 
at s = So (and might even have a zero there). 



Lemma 7. Let g{s) = L{s)f{s)j^. Vj G Z, set Oj = a + ^i. If 



1. Vj € Z \ {0}, / is analytic at s = Oj, 

2. 3A < 1, B > 0, such that for all integers positive integers q, \f^'^^{0j)\ = 
0(|_7|"^log^ (where the constant in the big O may depend upon q) 

3. \fj G Z, L(s) has a pole of order at most ui at s = Oj; 

furthermore, the coefficients of the Laurent series of L{s) are identical at 
each s = 9j, 

4- s(s+i) ^ P^^^ order at most n2 at s = 9q, 



then the sum of residues at s = 9j can be written in the form 

ni+n2 — 1 ni — 1 

5^i?es(5(s),s = e,)= A,n'^lg'n+ ^ Fi(lgnKlg^n, (21) 



«=ni 



i=0 



where the Aj 's are constants and Fi{u) 's are periodic functions with period one 
given by their Fourier series Fi{u) = Ylij^z'^id^'^^^^^ ■ Furthermore, all the 
Fourier series Fi{u) are absolutely convergent. 



Proof. We first introduce a notation. When r is clear from the context, 

[m I am, Qm+l, am+2, • • • ] 



represents the Laurent series Yli^o ^m+i{s — r 



\m+i 



We start by stating the Laurent series of each factor of g{s) at s 
ViGZ\{0}: 



0,-, where 



L{s) 
fis) 



s + 1 



-Ul 



lo, h,l2,h, ■ 



n = 

1 

s 



1 1 





1 1 



1 

1 1 
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The residue of g{s) at s = Oj is obtained by multiplying all these series together 
and extracting the coefficient of the term [s — Oj)^^- The residue will therefore 
be the sum of terms, each term of the form 



1 In^^ 2 
Ix, X —f(^^){ej) X n'^e^^^^'g'^ -Ig^ 



n X 



(-1) 



X4 



(-1) 



X5 



where Vi, Xj > and Ylf=i Xi = ni — 1. 

Hence the sum of these residues, when sorted according to the variable X3, 



IS 



ni — 1 



.2:3 



n] 



X3=0 

where 



ln^3 2 



E Jj(xi,X2,X4,X5,ni-l-X3)e2-^^''^" , (22) 



UGZ\{0} 



Jj(xi,X2,X4,X5,r) 



V il X J-f(-)f0 )xi^x 



+ 1)^5+1 



The relevant Laurent series of g{s) at s = 9q are: 



Lis) 



ni 


-n2 



lo, h,l2,h, 



1, Inn, — In^ n, — In^ n, • • • 
2! 3! 



ho,hi,h2,h3, 



s{s + l) 

By multiplying all these series together and extracting the coefficient of the 
term (s — Oq)~^, the residue at s = is found to be of the form 

ni+n2 — l rti — 1 

E Wlg*n+ E ^*^"lg'^- (23) 



i=ni 



i=0 



The second summation in (|23p combines with (|22p to give the second summation 
in ()2ip and the first summation in ()23p gives the first summation in (j2ip . 



We now prove the absolute convergence of the Fourier series. Take M = 
maxo<i<„i_i \li\. Note that for q<ni-l, \ f^'^\9j)\ = 0(|j|^log^ |j|). Thus 

|Jj(xi,X2,X4,X5,ni - 1 - X3)[ 
77-1 + 2 - X3 



< 



1 t ^ (-1)^'* 



< ('"^ + 3""^lxMxO(|j|^log^ 



xO( ^ 



O 



log 



|2-A 



Since {2 — A) > 1, the Fourier series is absolutely convergent. 



□ 
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To conclude, we note that as we only upper bound the order of poles but 
do not know their exact order, Aj may be zero and the Fi{u)^s may be constant 
functions, or even zero functions. 



3 Multidimensional Divide-and-Conquer 

3.1 Background of Multidimensional Divide-and-Conquer 

Multidimensional Divide-and-Conquer (MDC) was first introduced by Bentley 
and Shamos [5l H] in the context of solving multidimensional computational 
geometry problems. The generic idea is to solve a problem on n d-dimensional 
points by (i) first splitting the points into two almost equal subsets and solving 
the problem seperately on each subset, then (ii) taking all n points, projecting 
them down to (d — 1) dimensional space and solving the problem on the pro- 
jected set, and finally (iii) constructing a solution to the complete problem by 
intelligently combining the solutions to the 3 previously solved ones. The re- 
cursion bottoms out when the dimension d = 2, in which case a straightforward 
solution is given, or when n = 1, which has a trivial solution. 

The methodology can be applied to give good solutions for many problems, 
including the Empirical Cumulative Distribution Function (ECDF) problem, 
maxima, range searching, closest pair, and the all nearest neighbour problem. 

Of particular interest to us is the all-points ECDF problem in M.^ (ECDF- 
k). For two points x = {xi,X2,--- ,Xk), y = {yi,y2,--- ,Vk) G , we say x 
dominates y li Xi > yi for all 1 < i < /c. Given a set S" of n points in M*^, the 
rank of a point x is the number of points in S dominated by x. The ECDF-fc 
problem is to compute the rank of each point in S. 

When A; = 2, a slight modification of bottom-up mergesort will solve ECDF- 
2 in B(nlogn) time. Monier [18] proposed an MDC algorithm for solving 
ECDF-/c for larger A;, based on the description of Bentley [1]. Monier ana- 
lyzed the worst-case running time of this algorithm, T{n,k), described by the 
following recurrence: 



By translation into a combinatorial path-counting problem he derived the 
first order asymptotic of T{n, k). More specifically, he showed that, for fixed k, 



We will derive exact solutions for the ECDF-A; running time using Lemma 




(24) 



T{n,k) 



1 



nig'' ^n-h6(nlg^ ^n). 



{k-l)\ 
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[2] from [To]. To do so, we will have to slightly modify the case fc = 2 to have 
a more precise initial condition. In what follows we will denote T{n,k) by 
The recurrences corresponding to (|24p will be: 

" 1 0, n = 1 

where 

n — 1, k = 2. 



3.2 Deriving the DGF 

To use Lemma [2] first requires a better understanding of the DGF of A\/fll, 
which we denote by by Df^{s). Start by noting that, directly from the lemma, 



1 ~ AVe^ 
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One can work out directly that AVef = 1 while, for j > 2, AVe| = 0. Thus, 

1 ^AVe? 1 

3=1 3=1 



For /c > 3, 



' le^ = /2'-^ + l = AV/^-i + l, fori = l. 



Hence 



°° AVe^^ 



' ' AV/,^-1 + 1 + ^ .'''^^^ 



1-2-M ^ 7 

1 



1 - 2-** 1-2- 



Iterating the above recurrence with initial condition (|27p yields 
1 1 1 

^fk — ~x + T\ + 7l o-s\fc-l ■ 



1-2-'* (l-2-'')2 (1-2- 
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From Lemma [21 for A; > 1, 

fk _ ra r^+'°° 1 \ n'ds 

^" - 2vrzy3__ [2^^(l-2-^r) s{s + l) 



We note that Flajolet and Golin [10] explicitly solve the k = 2 boundary 
case: 

fn = n\gn + nAl{lgn) + 1 

where, setting f3j := ^^i, 



jez\{o} 



3.3 Evaluation of Integrals 



We now evaluate the integral in (|28p : 

3+joo 



n f^+"^ 1 n'ds 



27TiJs-ioo (l-2-«)'=-is(s + l)' 



(29) 



Fix some real R > and consider the counterclockwise rectangular contour 
T = Ti U T2 U T^s U ^^4, where (see Figure H 

Ti = {3 + iy: -R<y < R} T2 = {x + iR : -R < x < 3} 

T3 = {-R + iy : -R<y < R} T4 = {x - iR : -R < x < 3} (30) 

Denote the kernel of the integral in (I29p by Kk{s): 

Kkis) = 7 ri—r-. r. (31) 

^ (1 - 2-«)'=-is(s + 1) ^ ^ 

Note that 4 = limR^^oo 2^ /x, Kk{s)ds. 

We now show that, for g = 2, 3, 4, lim/j_>.oo /x Kk{s)ds = 0. Thus 

h = lim TT^ / Kk{s)ds. 

By the Cauchy residue theorem, Ij. will be equal to n times the sum of the 
residues at the poles inside T as -R — >• cx). 

The poles of K]^{s) inside T are: 
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Im{z) 
R 



Tj 



Ti 



Relz) 



Figure 4: Contour T defined in (j30|) . The dots represent the poles of Kk{s) 
inside T. 



1. A pole of order /c at s = 0; 

2. Poles of order (fc - 1) at s = (3j = ^i, where j G Z \ {0}; 

3. A simple pole at s = — 1. 

To avoid poles of Kk{s) on T, we only consider values oi R = Rj := ^'^^'1^2^^ ■ 
Now consider the horizontal paths q = 2,4. Then 



Kk{s)ds 



< 



< 



max 



\Kk{s)\ds 



(1 ± 2-~"i 



Rji^j + 1)/ J-R 



da = Oir') 



For the leftmost path it is easy to see 



Kk{s)ds 



O 



Hence Ik is n times the sum of the residues at the poles of Kk{s) inside T, 
taking Rj ^ oo. 

Theorem 2. For k > 2, 

k-2 

= 7^3W"lg'"'^+ E (^lg™n)^i(lgn) + Cfc, (32) 



m=0 
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where A'^{u) 's are periodic functions with period one, which are given by abso- 
lutely convergent Fourier series 

whose coefficients cik^mj can be determined explicitly. In particular, the average 
value of A'^ ^i.''^) 

1 fk + 1 1 



(A; - 2)! V 2 \n2j ' 
Furthermore, if k is even, Ck = 1; if k is odd, Ck = 0. 

Proof. We proceed by induction on k. As previously mentioned, for k = 2 this 
theorem was aheady proved by Flajolet and Gohn [lOj. 



Now assume that (j32|) is true for k = ko — 1. The residue of ii'fcp(s) at 
s = — 1 is 

(_l)fco 



Res{Kk^{s),s = -1) 



n 



We can now apply Lemma [71 by taking a = 0, L(s) = (1 — 2 (its 
Laurent series coefficients at each s = /3j are identical) and /(s) = 1. Since 
/(^^(s) = when q>l, we may take A = B = 0. The order of poles of L{s) at 
s = f3j is /co — 1 and the order of pole of J^^^^^^ at s = /3o = is 1. 

The sum of residues at s = /3j, where j G Z, is given by 

fco-2 

Afc„_ilg'=o-in+ (lg™n)S^o(lgn), 

m=0 

where i?^(u)'s are periodic functions with period one which are given by ab- 
solutely convergent Fourier series. Afc„_i can be explicitly calculated to be 
l/{ko - !)!• 



17 



Hence by (j28l). 

Jn 



I ft;n — I I : 

m=0 



+ j;(lg-nXo(lgn) + 



^ u ^ m=0 

+ ^(nlg'"n)<o-\lgn) + Cfc„_ 



nig " n 

1 AiQ — 1 ^o~2 

+ a^_lV +E (r^lg'"n)Sio(lgn) + (-l)^° 

° m=0 

feo-3 

+ E (^m"'(lg^) + <°(lg^)) nig- n + c,„_i + (-1)*^°. 

m=0 

Letting <_2(n) := + <'_2(^) and A^-^u) + B^iu) 

for m = 0, 1, • • • , fco — 2 proves (j32]) . The average value of A^_2(n) is found by 
expressing all the Laurent series (in the proof of Lemma [7]) explicitly. 

Finally, since = CfcQ-i + (— 1)'^", alternates between being even and odd 
with C2 = 1. □ 



4 Weighted Digital Sums of the First Type 



We now analyze TSM{n) = Ylj<n^Mij) as defined in ([2]) and ([3]). By Lemma 
[H this reduces to evaluating 

-MH^irte^);^. (33) 



C — lOO 



fr{ f i <s + 1) 



4.1 Deriving the DGF 

We start by deriving a closed form for 



Am{s) := 2^ . (34) 



Recall that Suin) = Y^\=Qt^^ bt2^ . Observe that if n = (6j?*j_i ••• 6160)2, 
then 

2n = (6i6i_i •• -6160,0)2 and 2n + 1 = (6j6i_i ••• 6160, 1)2. 
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In particular, when M > 1, the weight t^'^ for the rightmost digit {t = 0) is 
always zero, so 

5A/(2n + l) = 5Af(2n). (35) 

Next, observe that 

i i i 

Si{2n) = ^(t + l)6t2*+i =2^t6t2* + 2^6t2* = 25i(n) + 2n (36) 

t=0 t=0 t=0 

and for M > 2, 



5M(2n) = J^(t + l)^^5i2*+i 

i=0 



2 t^'h2' + M ^(t + l)^^-i6t2*+i 

t=0 t=Q 

2SM{n) + MSM-i{2n). (37) 



These facts lead to: 
Lemma 8. 



^Mis) = M!^,^^3^— ^C(s). (38) 



2{M-l){s-l) 

'(2^-1 - 1)^ 

Proof. The proof is by induction on M. When M = 1, by (j35p and (|36p . we get 
V5i(2n) = 2V5i(n) + 2 and V5i(2n + 1) = 0. 

Hence 

V^ZMZ) ^ 2V5i(/) + 2 1,,,,^., 

Then Ai[s) = [2^^^ — l)^^C(s) and the lemma is proved for M = 1. 
Now assume the lemma is true for M < k. Iterating (|37p gives 

5fc(2n) = 2Sk{n) + 2 k^Sk-i{n)^ + 2k\n, 

where k- = k{k — 1) ■ ■ ■ (k — i + 1) is the i^^ falling factorial of k. 
Appling ([35]) gives 

\/Sk{2n) = 2V5fc(n) + 2 ^ fc^V5fc_i(n) j + 2k\, 
VSk{2n + l) = 0. 
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Substituting the above two formulae into (jMj) yields 



Ms) = E 



and hence 



j=i 1=1 ^ ' 

2VSk{l) + 2 (Eti fe^V5fc_i(0 ) + 2A;! 



1=1 \ 



- 



{21Y 



1 ^ k \ ^1 



— Ms) + 



r^.(^) + l^(i + E^ 



2(fc-i-l)(s-l) 



E^-(fc 



« ! 



■t=l 



(2^-1 - 



C(s) + A:!C(s) 



fc-i 



2(A;~j-l)(s-l) 



— 1 ) ' 2^ \ 2^ ^ V 2'^ — 1 



k-l 



1=1 



2s-l^k[S)+ I 2^-1-1 



2(fc-i)(s-i) 



□ 



4.2 Evaluation of the Integral 

Substituting the result of Lemma [8] into the integral of (j33p gives 



Fix some real R > and consider the counterclockwise rectangular contour 
r = Ti U r2 U Ts U where (see Figure [5]) 

Fi = {3 + iy : -R < y < R}, T2 = {x + iR : -1/4 < x < 3}, 

Fg = {-1/4 + : -i? < 2/ < i?}, r4 = {x - ii? : -1/4 < x < 3} (40) 

Denote the kernel of the integral in (p9|) by KMis)- 

2{M-l){s-l)^( \ s 
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Im(z) 



r2 



1 + 

^ + ln2' 



r3' 



1-3,1 



3 



^ Re[z) 



1 - —7 



r4 



Figure 5: Contour F from (|40p . The dots represent the poles of Km{s) inside 
F. 

Note that TSM{n) = Hm/j^oo ^ /p^ KM{s)ds. As in the MDC case, we now 
show that liniR^oo /p KM{s)ds = for g = 2, 3, 4. Thus 



M! f 

TSM{n) = lim — / KM{s)ds. 



Hence by the Cauchy residue theorem, TSM{n) will be equal to Ml times the 
sum of the residues at the poles inside F as — > oo. 

We know that (^{s) has a simple pole at s = 1. The poles of Km{s) inside 
F are: 



1. A pole of order (M + 1) at s = 1; 

2. Poles of order M at s = := 1 + ^i, where j £Z\ {0}; 

3. A simple pole at s = 0. 



To avoid poles of Km{s) on F, we again only consider values oi R = Rj 

(2j+0.5)7r 
In 2 
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To show that /p KM{s)ds = for g = 2,3,4 as i? —> oo, we need the 
following two lemmas. 

Lemma 9. Consider integral 



a+iR 



f{s)as)n'ds, 



where < a < |. Furthermore, suppose that for s = a + it with —a < a < 3, 
|/(s)| = 0(1*1"^). Then, both as R ^ oo and R ^ -oo, I{R) 0. 



Proof. Along the path of the integral, 3?(s) > —5/4. Lemma [3] gives the bound 

\C{s)\=0{\Rf/Hog\R\). 

Together with the given fact that |/(s)[ = O 

r3+iR f3+iR 



-a+iR 



f{s)as)n'ds 



< 



\f{s)C{s)n'\ds 



i-S+iR 

< / (0{\R\-^) xO{\R\'^/^log\R\) 

J-a+iR. ^ 

< (3 + a)0(|i?r^/^log|/?|)n^ 




X n 



3^ ds 



as i? — 7- oo. 

Lemma 10. Suppose 



□ 



g{s) = f29j{K,y 

j=0 

for some real sequence {gj} and positive integer sequence {Kj}. If this series 
is uniformly convergent for sg{ — : t€M} then 

J-l/4.-ioo S{S + 1) 



If the series is uniformly convergent 

I — 5/4+joo 



/— 5/4— ioo 



for s e {-\ + it : t eM] then 
n^ds 



s{s + l){s + 2] 



= 0. 



Proof. For the first integral, note that 



-1/^+*°° , , n'ds 
9{s)C[s) 

— 1/4:— ioo 



s(s + l) 



-1/4+ioo / ~ » s J 







-^1^+'^ ^^^^ {Kjufds 

l/4-ioo ■s('S + 1) 
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The first equality is the definition of g{s), the second follows from the uniform 
convergence of the series and the last equality follows from (jl4p . 



The second integral is evaluated similarly, using p3|) in place of ([H]). □ 

To evaluate the integrals along Fa and note that |2(^-^)(^-i) (2*"^ - 1)"^^ 

0(j~^). Thus, by Lemma [H as Rj — > cxd, 



is bounded as j — )• oo and 



1 

s{s+l) 



To evaluate the integral along T^, note that c < along Fs, so we may 
write 



>s-l 



1 \2 V4 



The series is both absolutely convergent and uniformly convergent on —1/4+ 
(— (X),oo)i, so we may write (see [211 PP-74-75]) 

9(M-l)(s-l) °° 

for some {aj}, where this new series is again uniformly convergent on — 1/4 + 
{—oo,oo)i. Bv Lemma [TOl 



lim / KM{s)ds 0. 



We have successfully shown that the integrals along and F4 vanish as 

Rj ^ CO, and hence TSMin-) is Ml times the sum of the residues at the poles 
of Km{s) inside F, after taking Rj 00. 



Theorem 3. For M >l, 

TSuin) = -n\g^'n+ ^ (nlg-^n) FM,d(lgn) + (-1)^+^M!, (42) 

d=0 

where Fm4(v) 's are periodic functions with period one, which are given by ab- 
solutely convergent Fourier series 

whose coefficients fM,d,j can be determined explicitly. In particular, the average 
value of Fm,m-i{u) is 

M 

/m,m-i,o = [270 - 3 + (M - 2) In 2] « — - 0.915648M. 
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Proof. As shown, TSM{n) is M\ times the sum of residues at the poles oiKM{s) 
inside F as i? ^ oo. The residue of Km{s) at s = is 



Res(i^M(s),s = 0) = (-1) 



M+l 



By LemmaO we have the bound \C,{a + it)\ = 0{\tY log \t\) when a >1 — e 
for some sufficiently small e. By Lemma O IC'^'^H'^j)! — 0{\jY\og\j\) for any 
fixed positive integer q. 

In LemmaEl take cr = 1, L{s) = 2(M-i)(s-i)(2"-i - 1)"^ (its Laurent series 
coefficients at each s = aj are identical) and /(s) = Ci^)- From last paragraph 
we can take A = e and B = 1. The order of poles of L{s) at s = aj is M, and 
the order of pole of gj^j^i^ at s = oq = 1 is 1- 

The sum of residues at s = aj, where j G Z, is given by 

A/-1 



AAfnlg^^n+ Ig"' nj FAf,rf(lgn), 



d=0 



where i^M,d(^)'s are periodic functions with period one which are given by 
absolutely convergent Fourier series. Am can be explicitly calculated to be 
1/(2M!). 



Hence 



TSM{n) = Ml 



, M M-l 

nig n ^ . 



2M! 

M-l 



+ Y (nlg'^n) FM,d(}gn) + i-r 



M+l 



= ^nlg^^n+ Y (nlg'^ nj M\FM,dilgn) + {-1)^^+^M\. 

d=0 

Letting Fma{u) ■= M\FM,d{u) for d = 0, 1, • • • , M-l, proves (j42]). The average 
value of Fm,m-i{u) is found by expressing all the Laurent series (in the proof 
of Lemma [7]) explicitly. □ 



5 Weighted Digital Sums of the Second Type 

We now analyze TWi{n) = J2j<n ^i(i) ^ defined by (jl]) and ([5]). The analysis 
will be extended to TWM{n) for M > 1 in the next section. 

The general methodology used to analyze TWi{n) is the same as in the 
previous sections; use Lemma [1] to rewrite 
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The main difficulty tliat will be encountered is that the DGF here will not be 
"nice" enough to permit integrating the kernel directly. We will have to split 
the DGF into two parts, using the m = 1 case of ^ to evaluate the first part 
and the m = 2 case to evaluate the second part. 



5.1 Deriving the DGF 

Set 

to be the DGF of VWmO)- We start by deriving, for all M > 1, a formula for 
Bm{s) in terms of DGFs Vm{s) and Zm{s) introduced in Definition [21 We will 
then analyze the case M = 1 in this section, and leave the cases M > 2 to the 
next section. 

Lemma 11. 

2^ — 1 1 
Bm{s) = ^73^^m(s) - ^j— ^Zm(s). 

Proof. Observe that if n is expressed as n = 2*^ + 2*^ + ■ ■ ■ + 2**= with ii > i2 > 
• • • > ^fc > 0, then 



WMi2n) = ^t^^2^*+i = 2^t*^2^* = 2WM(n), 
t=i t=i 

k 

WM{2n + l) = Y.t^'2''+' + {k + lf' = 2WM{n) + {vin) + l) 



t=i 



Recalling from Lemma|3]that f(n)—t;(n—l) = 1—V2{n) andu(2n+l) = v{n) + l 
gives 



VWM{2n) = 2VWMin)-{v{n)+V2{n))^^ and VM^M(2n+l) = v{2n+l)^^ 

Then, (jlSp in Lemma [5] permits writing 

„ . ^ V- , 2VWM{l)-Hl)+V2{l)r 

Bm{s) = L-^ + L 

odd j -^1=1 ^ ^ 

Solving for Bm{s) proves the lemma. □ 
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For M = 1, applying (jl9p from Lemma [5] to Lemma [TT] gives 



Substituting this into (j43|) yields 

1 /-S+ioo s , 1 /.3+ioo A/ \ SI 



(46) 

The second integral can be evaluated exactly by the method used in Section 
Evaluating the first integral requires more work. 

Historically, v{n) was one of the first digital functions to be analyzed using 
the Mellin transform techniques. The original analysis in 1975 by Delange [8] 
used a combinatorial decomposition of the binary representations of integers to 
directly derive an exact Fourier series formula for Ylj<n'^U)- 1994, Flajolet 
et. al. [9] reproved Delange's result using the Mellin transform techniques. 
However, Vi{s), the DGF of v{n), does not seem to have been explicitly studied 
before Hwang's analysis [14] in 1998. First, denote Ir{s) by 

1 °° /I 9 1 \ 

/^.(s) := -iy 7;(i)M -V - ^ + 7 -^]- (47) 

^' V(22)^ (2^ + 1)^* (2i + 2)V ^ ^ 

Standard algebraic manipulations, e.g. in [HI pp.536], let us rewrite a summa- 
tion of this form as an integral: 

Iri^) = i f "-^a-)d-, (48) 

where 

-i ,if La;J <x< [x\ + i. 



e(x) 



,if [xj + i < X < [xj + 1. 



Hwang [14j derived the following formula of Vi{s), revealing its singularities 
in ^{s) > -1: 

V.is) = ^C(.) - ^(^C(.) + ^/i(.), (49) 
Substituting ([MD into (gS]) yields 
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From the integral form of Ii{s), we know that it is analytic in > — 1. 

Hence, (jSOp (together with the fact that (^{s) has no zero on the line K(s) = 0) 
shows that at s = 0, /3j, Bi{s) possesses simple poles. Depending upon the 
values of /i(l) and Ii{aj), Bi{s) may either possess simple poles at s = 1, aj 
or be analytic at s = 1, a^. These are all possible poles of Bi{s) inside T which 
we defined in (00]). 

Using Hwang's representation would yield a closed-form formula for TWi{n) 
by considering contour F. Unfortunately, the residues appearing in the resulting 
Fourier coefficients would be expressed in terms of the value of /i(s) at various 
poles, something which is not well understood. In the next subsection, we 
will show how to use the higher order version of the Mellin-Perron formula to 
sidestep this issue and express the Fourier coefficients in terms of the Riemann- 
Zeta function. 

5.2 Moving Up to a Higher Order Case of the Mellin-Perron 
Formula 

We now see how to manipulate the first integral in ()46p to yield a formula in 
terms of values of the Riemann-Zeta function. 

The general approach is to note that Vi{s) is the DGF of v{j), so the first 
integral in (j46p . when transformed from integral back to summation by (jlOp . 
is a double summation of v{j). A double summation of v{j) is also a triple 
summation of Vt'(j), and we can write a closed-form formula for the DGF of 
'Vv{j) in terms of C{s)- Equation ()lip then provides an exact formula of the 
triple summation of 'Vv{j), and we can evaluate the first integral in ()46p . 

We now present the details. Define 




n j-l 



j=l i=l 



Algebraic manipulations permit writing TV{n) in two different ways: 




(51) 



k<n 



and 



TV{n) = - ^S/v{k) 



(n — k)'^ + {n — k) 
2 



(52) 



k<n 



Applying ([TO]) to dH]), yields 




(53) 
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where the right side is exactly the first integral in (j46p . 

Applying (fTTj) and to (f52|) gives the alternate expression 

TV{n) (54) 

Setting ([53|) equal to and using the closed-form formula for ^v{3 )j^^ 

in Lemma [5] gives 

1 f^^'^yf n'ds 

n f3+ioo2s_2 n'ds 1 fS+ioo 2^ _ 2 n'ds 

-C(«)77:^^T7:^ + 7- / T^C(s)- 



27ri J3_,^ 2- - 1 " ^ ' s(s + + 2) 47ri ^g.^o, 2- - 1 " ^ ^ s(s + 1) 



Substituting the above equality into (j46|) yields a "nicer" integral represen- 
tation for TWi{n). 

TWiin) = — / — — C(s) 



27riJs^i^ 2-^-r' 's{s + l){s + 2) 
1 2^^ - 2 , . n'ds 



AiriJs^,^ 2'-V' 's{s + l) 
1 1 , n'ds 



C(^)^-rT- (55) 



27riy3-*oo (2^-l)(2^-2)^^ 's{s + l) 



5.3 Evaluation of Integrals 

The three integrals in (|55|) can be evaluated almost exactly as in Section 27 
That is, for the first integral consider contour T' = T[ U r2 U Fg U where 

r[ = {3 + iy: -R<y < R}, T'2 = {x + iR : -5/4 < x < 3}, 

r^j = {-5/4 + iy : -R<y < R}, = {x - ii? : -5/4 <x< 3}. (56) 

For the second and the third integrals consider contour T defined in (|4U|) . Next, 
prove that the integrals along the left, top and bottom paths tend to zero (using 
Lemma [9] and Lemma [TO]) . Finally, evaluate the sum of residues at the poles 
inside F' or F. Since these are almost exactly the same as in Section 14.21 we 
leave out the details, only stating the results. See Figure [6] for the contours. 

The poles of the kernel of the the first integral inside F' are a double pole at 
s = and simple poles at s = /3j (where j ^'L \ {0}) and s = —1. By summing 
the residues at all these poles, the first integral evaluates to 

n /■3+ioo2. _2 n'ds 1 n n 1 
/ C\s)—, r-, r = -nlgn + nili Ign H — , 57 
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Im{z) 



r'. 




. Re{z) 



Im{z) 

R r2 



I 1 
-II' • 1 - si' 

-eI' • 1 - ra' 

-i!l' • 1 - II' 



r4 



Figure 6: The coutours we use when evaluating the three integrals in (j55p . with 
poles marked with dots. The left one is T' defined in (|56p . which is for evaluating 
the first integral in (j55p . The middle and right ones are both T defined in (j40p . 
which are for evaluating the second and third integrals in ()55p respectively. 



where 
Hi{u) 



21n7r-ln2-3 \ 1 ^ C(/3j) „2..,. ..o^ 

81n2 ; ln2^.^^^^/3,(/3, + l)(/3, + 2) ' ^ > 



The poles of the kernel of the second integral inside T are a double pole at 
s = and simple poles at s = (3j (where j G Z \ {0}). By summing the residues 
at all these poles, the second integral evaluates to 

where 

/21n7r-hi2-2\ 1 ^ ({(3,) ,,,,, 
^^^^^"1 8hr2 J-2hr2 WTTiy" • 



The poles of the kernel of the third integral inside T are a double pole at 
s = 1, simple poles at s = aj (where j € Z \ {0}), a double pole at s = and 
simple poles at s = /3j (where j G Z \ {0}). By summing the residues at all 
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these poles, the third integral evaluates to 

f>3+joo 



V^h .00 (2- - 1)(2^ - 2) ^^^)^^ ^ ^nlgn+ng3,2(lgn) + -lgn+ff3,i(lgn 

(61) 



2'Ki 
where 



and 



/21n^+31n2-2\ 1 ^ C(/3,) o..,. 



' V 81n2 ; 2 In 2 ^.^A.^^ «,(«, + !) 



Combining the three integrals above yields: 
Theorem 4. 

TWi{n) = nFw,i{\gn) - ^ Ign + Fwfl{\gn). (64) 

where Fw,i{u) and Fwfi{u) are two absolutely convergent Fourier series, whose 
coefficients are given by 

ln7r-7o + 21n2 1 ^ ( 2C(/3j) C(aj) 
^'^^ ^ 41n2 21n2^.^Z^^^V/3,(/3, + l)(/3,+2) + «,(a,+l) 

_ 2-21n7r-51n2 1 ^ C(/?,) .a.^-. 
- 8h^2 + 2hr2 2- /3,(/3, + l)^ 

T/ie average value of Fw,i{u) is 

In vr - 70 + 2 In 2 
41n2 



0.704687. 



Proof. Substituting ([STD, ([59]) and dH]) into (l55|) yields 



TWi (n) = n [Hi (Ig n) - (Ig n)] - ^ Ig n + 



i/2(lgn) -i^3,i(lgri) + ^ 



Setting Fw,i{u) := Hi{u) - H^^2{u) and Fwfl{u) := H2{u) - Hs,i{u) + \ gives 
([Mj) and the Fourier series. 

Lemma [3] gives 

C(/3,) = 0(|ir/'log|i|) and C(a,) = 0(log |j|). 

Hence, as |j| 00 the terms in (|58l) . (i60]l and (|62]) are 0(|jj~^^^ log |j[) and the 
terms in ([S5|) are 0(|j|~^log implying the absolute convergences of Hi{u), 
H2{u), H?j^i{u) and H?j^2{u), and thus Fw,i{u) and Fwfi{u). □ 
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6 More Weighted Digital Sums of the Second Type 



We now analyze TW^f (n) = Y2j<n ^Af(j) as defined by ([6]) and d?]). Again, by 
Lemma [U 

1 fc+ioo sj 

where 5m (s) is the DGF of VWmU) as defined in (jM]). 

As before, we wih integrate along contour T we defined in (I40p . and prove 
that the integrals along the top, bottom and left contours vanish as R ^ oo, 
while that on the right contour equals ([65]) and then apply the Cauchy residue 
theorem. 

The DGF Bm{s) for Af > 1 is much more complicated than the DGFs 
previously encountered in this paper. We will therefore have to introduce new 
techniques to study it. 



6.1 Properties of Poles of the DGF 



We saw from (f50|) that the order of the poles of "^^^^^^"^ at s = 1 and s = aj 
were all at most 1. By Lemma [71 this implied that "coefficient" of the first 
order term of TWi{n), i.e. the n term in (j64p . is the Fourier series Fw^i{lgn). 

Analogously, for all M > 1 we will prove that the orders of poles of ^^^j^^ 
at s = 1 and s = aj are all at most 1. This will again imply that for M > 1, 
TWMin) has a "periodic first-order coefficient". 

To start, we will need the following semi-recursive formula of Bm{s)- 
Lemma 12. For M > 1, Bm{s) satisfies 



Bm{s) 



-as) + 



1 



1 ^ 

r=l 



2s ' 2 

where Rr{s) is defined as 

oo 

Rr{s) ■.= Y,V{ 



Br{s) 



r=l 



i=l 



1 



{2iY {2i + 1) 



M 



Rr{s), (66) 



(67) 



Proof. Lemma [TT] gives 

Bm{s) 



-Vm{s) 



-Zm{s). 



(68) 



We now derive two seperate functional equations for Va/('S) and Zm{s) and 
combine them to yield ()66p . 
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Recalling from Lemma U] that v{2i) = v{i) and v{2i + 1) = v{i) + 1 gives 
1 „ , > , /'"A /'■f^ "(i)' 1 



Solving for VMis) yields 

VMis, = l^^co) + ^ E ) v;(.) - ^ E (T)^'C)- 



Next, from Lemma[l]it is easy to show v{2i) + V2{2i) = v{i) + V2{i) + 1 and 
v{2i + 1) + V2{2i + 1) = v{2i + 1). Also, using (fTSjl in Lemma[5l gives 

^m(s) = 2^ T^T- + 2^ — ^ 

i>l ^ ^ odd j •' 



Since Zq[s) = C(s), this solves to 

Zm{s) = Vm{s) + ^C(^) + ^ E ( J 

r=l 



M-1 

/ ivi \ 

s). (70) 



Substituting ([69]) and ([70]) into ([68]) yields 

2^ 1 ^^7m\ (2--l)y,(.)-Z,(.) 2^ ^f^\j, 

r=l ^ ^ r=l ^ ^ 



(71) 

Finally, using (|68p to simplify the internal terms in ()7ip yields (166p . □ 



The next lemma gives a "closed- form" formula for Bm{s), in terms of /^(s), 
previously defined in (f^7|l . and C('5)- 



Lemma 13. 

M 



„ . X Pa/,i(2'') ^/ X , PM,2,k{2'') , . , . 

Bm{s) = _ ^^^, C{s) + (2^ _ i)M-fc(2s _ 2) ^^(^^' 

where Pm,i{x) and PM,2,k{x) cif^ iwo polynomials, with Pi^i(l) = i, Pm,i{^) = 
MPm-i,i{1) - M\/2^''^ for M>2 and PM,2,k{0) = for k = 1,2, ■■ ■ ,M. 
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Proof. We prove by induction, using Lemma [T2l First, note that (j50p is just 
the special case of this Lemma for M = 1. It also gives i-i^i(l) = ^. 



Now assume the lemma is true for M < Mq. Then 
^ Mo-i 

2 



r=l ^ ^ 

^ (2^ - i 2^* - 1 4^ (2'' - l)'-'=(2* - 2) ^ ^ 

I Z-,- (2'^ — J ' ^ -'fcwJ (2'^ — l)^+i~'^(2* — 2) 

._ rMo,i(2^) . . ^Y^^ ?^Afo,2,fc(2^) . . . . 

•~ (2'' - l)*^o ^ (2^ - l)^^o-'^(2'^ - 2) ^ 

where TA/f,,i(j;) and TMo,2,fc(a;)'s are polynomials satisfying rA/(,,i(l) = (jvi|f° J-Pa/q- 
MoPa/o-i,i(1) and rMo,2,fc(0) = for = 1, 2, • • • , Mq - L 

Now set 

1 ^ v[v{iY] 



From the definition of Ir(,s) in ()47p and some algebraic manipulation, 



Grabner and Hwang [13] proved 

2^^ ^ J^!5(rvfe^ i^/r\ A 

fc=l ^ h=\ ^ ^ fc=l ^ ' 

where S{n^ k) are the Stirling numbers of the second kind. Noting that S{n, n) 
1, Rr{s) can be rewritten as 

2^-2 0^X2^ ^^^2^ 



where Qr,i(a;) and (5r,2,A:(3;)'s are polynomials satisfying Q^^i^i) = '^^'^^r''^'^ 
r!/2''. 
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This permits writing 



2" - 2^ \ r 

r=l 



/Mo (A^o)Q^_^(2^)\ 2^ ^ ^ (*^")'3n2,fe(2') 

I (2'5 _ x)?- j 2^^ - 2 A^'^'"'-'^'' (V - l)'"-^ 

\ — 1 / 1 T" — k 

C/Mo,i(2') ... ^ , 2^ ^ ^Mo,2,fc(2^) . . . . 



(2^-l)*^o^^ ^ 2^ - 2 ^ (2^ - 

where ?7j\/,)^i(x) and f/Mo,2,fc(a^)'s are polynomials satisfying C/Afp^i(l) = Mo!/2*^''_ 
Substituting ([73]) and ([71|) back into ([66]) gives 

/ N 2^ / ?Afo,l(2'') ^ . 7A/o,2,fc(2^) r . A 
^Afis) = ^73T^(^) + I (2s_i)Mo ^(^) + (2^_i)Mo-fc(2^-2) ^^^^^ ) 

( UMoA'^n . , 2-[/Mo.2.fc(2'-) \ 

I (2^ - ^ (2^ - l)A'^o-fc(2s - 2)^^ ' j ' 

The lemma is proved for M = Mq by setting Pmo^i{x) := x(x — l)A^o~i + 
7Mo,i(^)-^A'/o,i(a^)> ^Mo,2,A:(a^) := TMo,2,fc(a;)-x[/Mo,2,fc(a;) for A; = 1,2,- •• ,Mo- 

1 and PAfo,2,Mo(a;) = -xUMo,2,Moix)- ^ 



We can now find the poles of Bm{s) inside T. See Figure [7] for locations. 

Corollary 1. For M >1, The singularities of Bm{s) inside F are 
(i) poles of order at most 1 at s = 1 and s = aj ; and 
(a) poles of order M at s = and s = (3j. 

Hence, the singularities of ^g^}^^ inside T are 
(i) poles of order at most 1 at s = 1 and s = aj; 
(a) a pole of order M + 1 at s = 0; and 
( Hi ) poles of order M at s = fij. 



Proof. (j72p permits us to identify the singularities by working through the 
various terms and recalling that Ik{s) is analytic when K(s) > —1. 

The recurrence relations Pm,i{^) = Af-PAf-i,i(l) — M!/2*^ with initial con- 
dition Pi,i(l) = 1/2 give Pm,i(1) > for M > 1. Hence at s = 0,/?^, 
-Pm,i(2*)/(2'^ — 1)*^ has poles of order exactly M, while C,{s) is analytic (but is 
not zero). 

At s = Uj, PAf,i(2'')/(2" - 1)^^ and C(s) are ah analytic. 
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, Re{z) 



Figure 7: The figure is contour F defined in (j40p . The dots represent the poles 

of ^M{s)r^ -^^-^^ r. 
s(s+l} 



At s = 1, Pm,i(2'*)/(2'* — 1) is analytic, but C{s) has a simple pole. 

At s = 0,Pj, the order of poles of PAf,2,fc(2")(2^ - l)-(A^-fc)(2^ - 2)-^Ik(s) 
is at most M — 1. 



At s = l,Qj, PM,2,fc(2^)(2^ - l)-(*^-'=)(2^ - 2)-i4(s) has poles of order at 



most 1 (due to the term (2* — 2) 



□ 



6.2 A Formula for rW^Af(n) 



As in the previous problems, we must again first show that the integrals along 
the top, bottom and left contours vanish as i? — t- oo. 

We need two basic observations. Suppose H{s) = P{2'^){2'^ — l)~^i(2* — 
2)~^'^, where P is a polynomial and Ni,N2 are non-negative integers. 

Fact 1: When 5R(s) < 0, H{s) can be expressed as a power series of 2^*, and this 
series is absolutely and uniformly convergent on the line 9?(s) + (— oo, +oo)i. 
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Furthermore, if P{0) = 0, i.e. the constant term of P is zero, then the constant 
term of the power series is also zero. 

Fact 2: H(s) is bounded along the line segment (—1/4, 3) + iRj independently 
of j. 



Lemma 14. 



j^ooJ_i/4+j^^ S(S + 1) 



Proof. For s S (—1/4,3) +iRj, Grabner and Hwang [13] proved that 

\lM{s)\ = 0{\jf /Hog''' \j\) = 0{\j\). 

Furthermore, Lemma [3] gives 

\Cis)\ = Oi\jf/Hog\j\)=oi\j\). 



By dZl]) and Fact 2, 1Bm(s)1 is bounded by o{\j\) along (-1/4,3) +iRj. 
Hence 



S(S + 1) 



< 



< 







s(s + l) 
(o(Iil) xO(ljr') xn^) 



as J — )• oo. 

Lemma 15. For any positive integer M , 



Bm{s)^—- = 0. 

l/4-ioo ■SI'S + -Lj 



□ 



Proof. Grabner and Hwang |13) proved the bound 



0(|^|3/4log2-|^|) = 0(|^|3/4+<5^ 



for any 5 > 0. This upper bound allows us to use a theorem from Hwang 
to prove 

Hwang [14] proved the following theorem: 

Theorem 5. Suppose U{s) = s2~* u{x)^{x)x~'^~^dx for some nonnegative, 
real arithmetic function u{x) = u^^.] ■ U 
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1. U{s) converges for K(s) > cr„, where au < cr, 

2. \U{a + it)\ = 0{\tf) for some < 5 < 1, 



then we have 

2TTi 



Ja-ioc S{S + 1) 



for all integers k,n>l. 

Grabner and Hwang [T3] proved the bound 

Ir ("i +^*) = Oi\tf /Hog'^ \t\) = 0{\tf/'+') 

for any 6 > 0, which enables us to use Theorem [5] to get 

±—l-Uis)ds = (75) 

for positive integers k,n,r. 

(j72]) shows that Bm{s) can be expressed in the form of 

^ ^ _ Pm,i(2-) , , ^ PM,2,fc(2-) 
^Ansj - (-2. _ i^M^y^) + 2^ (2^ _ i)M-fc(2« - 2) 

while PM,2,fc(0) = 0. By Fact 1, when 3?(s) = -1/4, Pm,i(2^)(2" - 1)-^^ and 
-PM,2,fc(2*)(2'* - 1)"(^^"'=)(2* — 2)^^ can be expressed as power series of 2^, and 
the power series for Pjvf,2,A:(2'^)(2'^ — 1)^'-*^~^^(2'' — 2)"-*^ have zero constant terms. 
Hence, when = —1/4, we may rewrite Bm{s) to be 

oo M oo 

bm{s) = Y.pj'^'Xis) + E E <ik,j'^''ik{s) 

j=0 k=l j=l 

for some {pj} and {qk,j}- Hence 

Bm{s) 



l/4-joo ■^(■S + 1) 

/■-1/4+ioo / oo M oo \ 

/ E^'^(2^'^)^c(«) + EE9fc.^(2^'™)%(.) 



l/4-ioo ^^^.^^ ' I s{s + l) 

However, the power series X]jlo^'j(^"'"')'* ^^^'^ Z^j^i ^fejC^"'"-)* uniformly 
convergent on —1/4 + (— cxD,oo)i, by Fact 1. This allows interchange of the 
integral sign and the summation signs. 

Hence, f_il^^-^ Bm{s) can be expressed as a series, in which each 

term is either a constant multiplied by an integral in the form of (|75p . or a 
constant multiplied by an integral in the form of (|14p . □ 
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We can now state our final result. 
Theorem 6. 

M-l 

TWM{n) = nGA/(lgn) + dM Ig*' n + ^ (ig-^n) GA/,d(lgn), (76) 

where dM is a constant, Gm{u) and GM,d{u) 's are periodic functions with period 
one given by absolutely convergent Fourier series. 



Proof. Consider the contour F in Figure [71 taking R ^ oo. Lemma [HI and 

r BM(s)n^ 

)Vq S{S+1) 



Lemma [T5] show that 2^ Jp ^s{s+t) '^^ = for (7 = 2, 3, 4. Hence 



TWM{n) = — / BM{s)^—rds, 
2m Jy^ s{s + 1) 



is the sum of residues at the poles of ^^^^^^^^ inside F, by the Cauchy residue 
theorem. 

By Lemma[3l we have the bound |C((T+it)| = Odtl"*^/^^*^ log \t\) when a > —e 
for sufficiently small e. Grabner and Hwang |13j also proved that \Ir{<7 + it)\ = 
0(|t|^/^''''^ log^*" when a > —e for sufficiently small e. Hence by Lemma[6l 

|C(^)(a,)l, lC('^^(/3,)l = 0(|jf/'+^log|j|) 

and 

|4^)(/3,)l=0(|jf/2+Mog2nj|) 
for any fixed positive integer q. 

Bm{s) can be expressed in the form of (|72p . Knowing that each function in 
the form of P(2*)(2* — l)~^i(2'* — 2)"^^ will have a Laurent series with identical 
coefficients at 6j = a + j^i for any fixed a, togather with the results from the 
last paragraph and Corollary [TJ we use Lemma [3 when cr = 0, 1 to derive 

VRcs ( ^^r^7^,s = aj] = nGAf(lgn) 

and 

^ Res ^f^, s = f3A= dM Ig^'n+Y^ (ig^ n) GmA^S n), 

where GMiu) and GM,diu)^s are periodic functions with period one given by 
absolutely convergent Fourier series. □ 
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7 Conclusion 



Mellin Transform techniques have previously been extensively used to analyze 
various divide-and-conquer algorithms and digital sums. A common theme in 
those analyses is the appearance of a (usually second order) periodic term, 
usually expressed in terms of a Fourier series. This Fourier series is the sum of 
residues of a complex function which has singularities regularly spaced along a 
vertical line. 

In this paper we pushed the technique further to derive exact analyses of 
the solution to multidimensional divide-and-conquer recurrences and various, 
more complicated, weighted digital sums. Our closed form solutions had the 
properties that a// terms were either polylogarithmic or n times a polylogarithm, 
with all coefficients either being constant or a periodic function given by a 
Fourier series. 

Our analysis of the multidimensional divide-and-conquer recurrence was a 
straightforward extension of the use of Mellin transform techniques for the 
analysis of simple divide-and-conquer recurrences. Our analyses of weighted 
digital sums, though, required developing a better understanding of various 
Dirichlet generating functions of differences of digital functions. 
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A Proof of Lemma [2] 

Our main technique for solving the Multidimensional Divide-and-Conquer re- 
currence is a generalization of Lemma [2] for basic divide-and-conquer recur- 
rences, originally proved in [TO] by Flajolet and Golin. In order to make this 
paper self-contained, we provide the proof of that lemma here. 

The divide-and-conquer recurrence is 

fn = /[n/2j + f\n/2\ + 

with initial condition /i = and given "conquer" cost sequence {e„} where 
eg = ei = 0. 

Distinguishing between odd and even cases of the recurrence, we find that 
for j > 1, 

fij = 2/j- -I- e2j, /2j+i = fj + /i+i + (77) 

Let V(7n = Qn — Qn-i be the backward difference operator. Then, for j > 1, 
V/2i = Vfj + Ve2„ V/2,+1 = V/,+1 + Ve2j+i. (78) 

Let /S.gn = Qn+i — dn, be the forward difference operator, i.e., 



(79) 



fn — V/„,+ i — V/„ — fn+1 — 2fn + fn-1 

AVe„ = Ve„+i - Ve„ = e„+i - 2e„ + e„_i. 

Then, from (j78p . 

AV/2J = AVfj + AVeaj, AV/a^+i = AVe2j+i 
for j > 1, with AV/i = /2 - 2/i = 62 = AVei. 

Basic calculations now show that, for any sequence /„, 

fn-nf, = ^{n-j)AVfj. (80) 

j<n 

Therefore, (llOp gives that 

fc+ioo SJ 
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where Df{s) := AV/jj'' is the DGF of AVfj. 

Further calculation yields 



00 



L>/(s) ^ AVej 



7 



Solving for Df{s) gives 



^ ^ i=i 

Combining (|81|) and (|82|) proves Lemma [21 

B Proof of dlS]) 

In this section, we mimic the proof of (jl4p in [9j to prove (|15p . 
Setting Xj = 1 in (fTT]l gives 

(n-l)(2n-l) 1 f^^'°° . n'ds 



, C{s)—, TT r- (83) 

Now consider the rectangular contour T', which we defined in (j56|) (see 
Figure ED. By ([MI), 

1 /■ . n'^ds (n - l)(2n - 1) 



lim / C,{s 



R^oo2m Jr'^ s(s + l)(s + 2) 12n 



By Lemma [9l 



C.{s) , , , and — / C{s)- 



2vri ir'a s(s + l)(s + 2) 27ri ^s(s + l)(s + 2) 

vanish as i? — > 00. The poles and their residues inside the contour can be easily 
computed. They are 
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I. 



-1 n 



7 



Figure 8: The contour V used to derive (fTSl) with the poles of g(g_(_i)(s_[_2) 
r' at -1,0, 1 noted. 



1. A simple pole at s = 1. The residue at this pole is ^. 

2. A simple pole at s = 0. The residue at this pole is ^C(O) = ~\ 

3. A simple pole at s = — 1. The residue at this pole is — ^tiLi = 



1 

I2n- 



By the Cauchy residue theorem, we have 



I — 5/4+joo 



27ri 



5/4— ioo 



n^ds 



lim — : / ({s) 



n 

6 + 



+ l)(s + 2) 
s(s + l)(s + 2) 



n 
6 
0. 



+ 



1\ ^ 
'4 j 12n 



1\ ^ 
4 j 12n 



y lim — /" C(i 

=1,2,4 



n^ds 



s{s + l){s + 2) 



(n - l)(2n - 1) 
12^1 
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C Proof of Lemma [5] 



Recall from Lemma U] that v{2n) = v{n), then 



odd j 



This yields 



js ^ i 

odd j i=i 

odd j 



odd j •' ^ ^ 



Recall from Lemma|4]that V2{2n) = V2{n) + 1; also, for all odd n, V2{n) = 0. 
We have 

V2{j) ^ ^'2(2i) 



2" ^ i 

1=1 



1 ^ V2ii) 



2s IS 



and hence 



Finally recalling from Lemma H] that Vf(n) = v(n) — v{n — 1) = 1 — V2{n) 
yields 
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